Deep resequencing identifies candidate functional genes in leprosy GWAS loci.

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.

Human genetics of mycobacterial disease.

Mycobacterial diseases are caused by members of the genus Mycobacterium, acid-fast bacteria characterized by the presence of mycolic acids within their cell walls. Claiming almost 2 million lives every year, tuberculosis (TB) is the most common mycobacterial disease and is caused by infection with M. tuberculosis and, in rare cases, by M. bovis or M. africanum. The second and third most common mycobacterial diseases are leprosy and buruli ulcer (BU), respectively. Both diseases affect the skin and can lead to permanent sequelae and deformities. Leprosy is caused by the uncultivable M. leprae while the etiological agent of BU is the environmental bacterium M. ulcerans. After exposure to these mycobacterial species, a majority of individuals will not progress to clinical disease and, among those who do, inter-individual variability in disease manifestation and outcome can be observed. Susceptibility to mycobacterial diseases carries a human genetic component and intense efforts have been applied over the past decades to decipher the exact nature of the genetic factors controlling disease susceptibility. While for BU this search was mostly conducted on the basis of candidate genes association studies, genome-wide approaches have been widely applied for TB and leprosy. In this review, we summarize some of the findings achieved by genome-wide linkage, association and transcriptome analyses in TB disease and leprosy and the recent genetic findings for BU susceptibility.

Genome-wide analysis of the genetic regulation of gene expression in human neutrophils.

Neutrophils are an abundant immune cell type involved in both antimicrobial defence and autoimmunity. The regulation of their gene expression, however, is still largely unknown. Here we report an eQTL study on isolated neutrophils from 114 healthy individuals of Chinese ethnicity, identifying 21,210 eQTLs on 832 unique genes. Unsupervised clustering analysis of these eQTLs confirms their role in inflammatory responses and immunological diseases but also indicates strong involvement in dermatological pathologies. One of the strongest eQTL identified (rs2058660) is also the tagSNP of a linkage block reported to affect leprosy and Crohn's disease in opposite directions. In a functional study, we can link the C allele with low expression of the ß-chain of IL18-receptor (IL18RAP). In neutrophils, this results in a reduced responsiveness to IL-18, detected both on the RNA and protein level. Thus, the polymorphic regulation of human neutrophils can impact beneficial as well as pathological inflammatory responses.

Understanding inflammatory bowel disease via immunogenetics.

The major inflammatory bowel diseases, Crohn's disease and ulcerative colitis, are both debilitating disorders of the gastrointestinal tract, characterized by a dysregulated immune response to unknown environmental triggers. Both disorders have an important and overlapping genetic component, and much progress has been made in the last 20 years at elucidating some of the specific factors contributing to disease pathogenesis. Here we review our growing understanding of the immunogenetics of inflammatory bowel disease, from the twin studies that first implicated a role for the genome in disease susceptibility to the latest genome-wide association studies that have identified hundreds of associated loci. We consider the insight this offers into the biological mechanisms of the inflammatory bowel diseases, such as autophagy, barrier defence and T-cell differentiation signalling. We reflect on these findings in the context of other immune-related disorders, both common and rare. These observations include links both obvious, such as to pediatric colitis, and more surprising, such as to leprosy. As a changing picture of the underlying genetic architecture emerges, we turn to future directions for the study of complex human diseases such as these, including the use of next generation sequencing technologies for the identification of rarer risk alleles, and potential approaches for narrowing down associated loci to casual variants. We consider the implications of this work for translation into clinical practice, for example via early therapeutic hypotheses arising from our improved understanding of the biology of inflammatory bowel disease. Finally, we present potential opportunities to better understand environmental risk factors, such as the human microbiota in the context of immunogenetics.

Association of a mutation in LACC1 with a monogenic form of systemic juvenile idiopathic arthritis.

OBJECTIVE: The pathologic basis of systemic juvenile idiopathic arthritis (JIA) is a subject of some controversy, with evidence for both autoimmune and autoinflammatory etiologies. Several monogenic autoinflammatory disorders have been described, but thus far, systemic JIA has only been attributed to a mutation of MEFV in rare cases and has been weakly associated with the HLA class II locus. This study was undertaken to identify the cause of an autosomal-recessive form of systemic JIA. METHODS: We studied 13 patients with systemic JIA from 5 consanguineous families, all from the southern region of Saudi Arabia. We used linkage analysis, homozygosity mapping, and whole-exome sequencing to identify the disease-associated gene and mutation. RESULTS: Linkage analysis localized systemic JIA to a region on chromosome 13 with a maximum logarithm of odds score of 11.33, representing the strongest linkage identified to date for this disorder. Homozygosity mapping reduced the critical interval to a 1.02-Mb region defined proximally by rs9533338 and distally by rs9595049. Whole-exome sequencing identified a homoallelic missense mutation in LACC1, which encodes the enzyme laccase (multicopper oxidoreductase) domain-containing 1. The mutation was confirmed by Sanger sequencing and segregated with disease in all 5 families based on an autosomal-recessive pattern of inheritance and complete penetrance. CONCLUSION: Our findings provide strong genetic evidence of an association of a mutation in LACC1 with systemic JIA in the families studied. Association of LACC1 with Crohn's disease and leprosy has been reported and justifies investigation of its role in autoinflammatory disorders.

Genetics of leprosy: expected and unexpected developments and perspectives.

A solid body of evidence produced over decades of intense research supports the hypothesis that leprosy phenotypes are largely dependent on the genetic characteristics of the host. The early evidence of a major gene effect controlling susceptibility to leprosy came from studies of familial aggregation, twins, and Complex Segregation Analysis. Later, linkage and association analysis, first applied to the investigation of candidate genes and chromosomal regions and more recently, to genome-wide scans, have revealed several leukocyte antigen complex and nonleukocyte antigen complex gene variants as risk factors for leprosy phenotypes such as disease per se, its clinical forms and leprosy reactions. In addition, powerful, hypothesis-free strategies such as Genome-Wide Association Studies have led to an exciting, unexpected development: Leprosy susceptibility genes seem to be shared with Crohn's and Parkinson's diseases. Today, a major challenge is to find the exact variants causing the biological effect underlying the genetic associations. New technologies, such as Next Generation Sequencing that allows, for the first time, the cost and time-effective sequencing of a complete human genome, hold the promise to reveal such variants. Strategies can be developed to study the functional effect of these variants in the context of infection, hopefully leading to the development of new targets for leprosy treatment and prevention.

Genetics of leprosy: expected and unespected developments and perspectives

A solid body of evidence produced over decades of intense research supports the hypothesis that leprosy phenotypes are largely dependent on the genetic characteristics of the host. The early evidence of a major gene effect controlling susceptibility to leprosy came from studies of familial aggregation, twins, and Complex Segregation Analysis. Later, linkage and association analysis, first applied to the investigation of candidate genes and chromosomal regions and more recently, to genome-wide scans, have revealed several leukocyte antigen complex and nonleukocyte antigen complex gene variants as risk factors for leprosy phenotypes such as disease per se, its clinical forms and leprosy reactions. In addition, powerful, hypothesis-free strategies such as Genome-Wide Association Studies have led to an exciting, unexpected development: Leprosy susceptibility genes seem to be shared with Crohn's and Parkinson's diseases. Today, a major challenge is to find the exact variants causing the biological effect underlying the genetic associations. New technologies, such as Next Generation Sequencing that allows, for the first time, the cost and time-effective sequencing of a complete human genome, hold the promise to reveal such variants. Strategies can be developed to study the functional effect of these variants in the context of infection, hopefully leading to the development of new targets for leprosy treatment and prevention.

M. paratuberculosis and Parkinson's disease--is this a trigger.

Genetic linkage studies and genome wide analysis have provided insights into complex medical diseases. Mycobacterium avium ss. paratuberculosis (MAP) causes Johne's disease, an important enteric inflammatory disease mostly studied in ruminant animals. MAP is also the putative cause of Crohn's disease. Moreover, MAP has been linked to other inflammatory diseases: sarcoidosis, Blau syndrome, autoimmune diabetes, autoimmune thyroiditis and multiple sclerosis. Genetic studies reveal an association between Parkinson's disease (PD), leprosy and Crohn's disease and since discovered, these findings have been considered "surprising". Autophagy and ubiquitin-proteosome systems are cellular systems that both fight intracellular pathogens (xenophagy) and maintain cellular protein quality control. PD is a common neurodegenerative disease that manifests clinically as a profound movement disorder. The recognized genetic defects of PD create disruption of cellular homeostasis that result in protein folding abnormalities of PD called Lewy bodies. Those same genetic defects are associated with susceptibility to intracellular pathogens, including mycobacteria. It is now understood that PD Lewy body pathology starts in the enteric nervous system and "spreads" to the brain in a retrograde fashion via the vagus nerve. This is the same process by which prions affect the brain. Lewy body pathology of the enteric nervous system predates the Lewy body pathology of the central nervous system (CNS) by years or even decades. This article proposes that genetic defects associated with PD also result in a permissive environment for MAP infection--ineffective xenophagy. It postulates that beginning as an enteric infection, MAP--via the vagus nerve--initiates a pathologic process that results in a targeted neuroinvasion of the CNS. The article proposes that MAP infection and resultant PD pathology are due, in the genetically at-risk and age dependant, to the consumptive exhaustion of the protein quality control systems.

A general efficient and flexible approach for genome-wide association analyses of imputed genotypes in family-based designs.

Genotype imputation is a critical technique for following up genome-wide association studies. Efficient methods are available for dealing with the probabilistic nature of imputed single nucleotide polymorphisms (SNPs) in population-based designs, but not for family-based studies. We have developed a new analytical approach (FBATdosage), using imputed allele dosage in the general framework of family-based association tests to bridge this gap. Simulation studies showed that FBATdosage yielded highly consistent type I error rates, whatever the level of genotype uncertainty, and a much higher power than the best-guess genotype approach. FBATdosage allows fast linkage and association testing of several million of imputed variants with binary or quantitative phenotypes in nuclear families of arbitrary size with arbitrary missing data for the parents. The application of this approach to a family-based association study of leprosy susceptibility successfully refined the association signal at two candidate loci, C1orf141-IL23R on chromosome 1 and RAB32-C6orf103 on chromosome 6.

CUBN and NEBL common variants in the chromosome 10p13 linkage region are associated with multibacillary leprosy in Vietnam.

Leprosy is caused by infection with Mycobacterium leprae and is classified clinically into paucibacillary (PB) or multibacillary (MB) subtypes based on the number of skin lesions and the bacillary index detected in skin smears. We previously identified a major PB susceptibility locus on chromosome region 10p13 in Vietnamese families by linkage analysis. In the current study, we conducted high-density association mapping of the 9.5 Mb linkage peak on chromosome region 10p13 covering 39 genes. Using leprosy per se and leprosy subtypes as phenotypes, we employed 294 nuclear families (303 leprosy cases, 63 % MB, 37 % PB) as a discovery sample and 192 nuclear families (192 cases, 55 % MB, 45 % PB) as a replication sample. Replicated significant association signals were revealed in the genes for cubilin (CUBN) and nebulette (NEBL). In the combined sample, the C allele (frequency 0.26) at CUBN SNP rs10904831 showed association [p = 1 × 10(-5); OR 0.52 (0.38-0.7)] with MB leprosy only. Likewise, allele T (frequency 0.42) at NEBL SNP rs11012461 showed association [p = 4.2 × 10(-5); OR 2.51 (1.6-4)] with MB leprosy only. These associations remained valid for the CUBN signal when taking into account the effective number of tests performed (type I error significance threshold = 2.4 × 10(-5)). We used the results of our analyses to propose a new model for the genetic control of polarization of clinical leprosy.

Chromosome 2p14 is linked to susceptibility to leprosy.

BACKGROUND: A genetic component to the etiology of leprosy is well recognized but the mechanism of inheritance and the genes involved are yet to be fully established. METHODOLOGY: A genome-wide single nucleotide polymorphism (SNP) based linkage analysis was carried out using 23 pedigrees, each with 3 to 7 family members affected by leprosy. Multipoint parametric and non-parametric linkage analyses were performed using MERLIN 1.1.1. PRINCIPAL FINDINGS: Genome-wide significant evidence for linkage was identified on chromosome 2p14 with a heterogeneity logarithm of odds (HLOD) score of 3.51 (rs1106577) under a recessive model of inheritance, while suggestive evidence was identified on chr.4q22 (HLOD 2.92, rs1349350, dominant model), chr. 8q24 (HLOD 2.74, rs1618523, recessive model) and chr.16q24 (HLOD 1.93, rs276990 dominant model). Our study also provided moderate evidence for a linkage locus on chromosome 6q24-26 by non-parametric linkage analysis (rs6570858, LOD 1.54, p = 0.004), overlapping a previously reported linkage region on chromosome 6q25-26. CONCLUSION: A genome-wide linkage analysis has identified a new linkage locus on chromosome 2p14 for leprosy in Pedigrees from China.

The Maximum-Likelihood-Binomial method revisited: a robust approach for model-free linkage analysis of quantitative traits in large sibships.

Model-free linkage analysis methods, based on identity-by-descent allele sharing, are commonly used for complex trait analysis. The Maximum-Likelihood-Binomial (MLB) approach, which is based on the hypothesis that parental alleles are binomially distributed among affected sibs, is particularly popular. An extension of this method to quantitative traits (QT) has been proposed (MLB-QTL), based on the introduction of a latent binary variable capturing information about the linkage between the QT and the marker. Interestingly, the MLB-QTL method does not require the decomposition of sibships into constituent sibpairs and requires no prior assumption about the distribution of the QT. We propose a new formulation of the MLB method for quantitative traits (nMLB-QTL) that explicitly takes advantage of the independence of paternal and maternal allele transmission under the null hypothesis of no linkage. Simulation studies under H0 showed that the nMLB-QTL method generated very consistent type I errors. Furthermore, simulations under the alternative hypothesis showed that the nMLB-QTL method was slightly, but systematically more powerful than the MLB-QTL method, whatever the genetic model, residual correlation, ascertainment strategy and sibship size considered. Finally, the power of the nMLB-QTL method is illustrated by a chromosome-wide linkage scan for a quantitative endophenotype of leprosy infection. Overall, the nMLB-QTL method is a robust, powerful, and flexible approach for detecting linkage with quantitative phenotypes, particularly in studies of non Gaussian phenotypes in large sibships.

Leprosy and the human genome.

Despite the availability of effective treatment for several decades, leprosy remains an important medical problem in many regions of the world. Infection with Mycobacterium leprae can produce paucibacillary disease, characterized by well-formed granulomas and a Th1 T-cell response, or multibacillary disease, characterized by poorly organized cellular infiltrates and Th2 cytokines. These diametric immune responses confer states of relative resistance or susceptibility to leprosy, respectively, and have well-defined clinical manifestations. As a result, leprosy provides a unique opportunity to dissect the genetic basis of human in vivo immunity. A series of studies over the past 40 years suggests that host genes influence the risk of leprosy acquisition and the predilection for different clinical forms of the disease. However, a comprehensive, cellular, and molecular view of the genes and variants involved is still being assembled. In this article, we review several decades of human genetic studies of leprosy, including a number of recent investigations. We emphasize genetic analyses that are validated by the replication of the same phenotype in independent studies or supported by functional experiments demonstrating biological mechanisms of action for specific polymorphisms. Identifying and functionally exploring the genetic and immunological factors that underlie human susceptibility to leprosy have yielded important insights into M. leprae pathogenesis and are likely to advance our understanding of the immune response to other pathogenic mycobacteria. This knowledge may inform new treatment or vaccine strategies for leprosy or tuberculosis.

Hereditary sensory neuropathy type I.

Hereditary sensory neuropathy type I (HSN I) is a slowly progressive neurological disorder characterised by prominent predominantly distal sensory loss, autonomic disturbances, autosomal dominant inheritance, and juvenile or adulthood disease onset. The exact prevalence is unknown, but is estimated as very low. Disease onset varies between the 2nd and 5th decade of life. The main clinical feature of HSN I is the reduction of sensation sense mainly distributed to the distal parts of the upper and lower limbs. Variable distal muscle weakness and wasting, and chronic skin ulcers are characteristic. Autonomic features (usually sweating disturbances) are invariably observed. Serious and common complications are spontaneous fractures, osteomyelitis and necrosis, as well as neuropathic arthropathy which may even necessitate amputations. Some patients suffer from severe pain attacks. Hypacusis or deafness, or cough and gastrooesophageal reflux have been observed in rare cases. HSN I is a genetically heterogenous condition with three loci and mutations in two genes (SPTLC1 and RAB7) identified so far. Diagnosis is based on the clinical observation and is supported by a family history. Nerve conduction studies confirm a sensory and motor neuropathy predominantly affecting the lower limbs. Radiological studies, including magnetic resonance imaging, are useful when bone infections or necrosis are suspected. Definitive diagnosis is based on the detection of mutations by direct sequencing of the SPTLC1 and RAB7 genes. Correct clinical assessment and genetic confirmation of the diagnosis are important for appropriate genetic counselling and prognosis. Differential diagnosis includes the other hereditary sensory and autonomic neuropathies (HSAN), especially HSAN II, as well as diabetic foot syndrome, alcoholic neuropathy, neuropathies caused by other neurotoxins/drugs, immune mediated neuropathy, amyloidosis, spinal cord diseases, tabes dorsalis, lepra neuropathy, or decaying skin tumours like amelanotic melanoma. Management of HSN I follows the guidelines given for diabetic foot care (removal of pressure to the ulcer and eradication of infection, followed by the use of specific protective footwear) and starts with early and accurate counselling of patients about risk factors for developing foot ulcerations. The disorder is slowly progressive and does not influence life expectancy but is often severely disabling after a long duration of the disease.

Genomewide linkage analysis of the granulomatous mitsuda reaction implicates chromosomal regions 2q35 and 17q21.

The Mitsuda reaction, a delayed granulomatous skin reaction elicited by the intradermal injection of heat-killed Mycobacterium leprae, is an in vivo test reflecting the ability to generate an immune granuloma after sensitization by diverse mycobacterial infections. Accumulating evidence for the genetic control of the Mitsuda reaction has been reported. We performed a genomewide linkage scan for the quantitative Mitsuda reaction in 19 large families from Vietnam with a history of leprosy (114 offspring). Suggestive linkage was found at chromosomal regions 2q35 (P = 9 x 10(-4) at the SLC11A1 locus) and 17q21-25 (P = 8 x 10(-4)). Interestingly, these 2 regions have been previously linked to mycobacterial infection and other granulomatous diseases.

Quantifying genomic imprinting in the presence of linkage.

Genomic imprinting decreases the power of classical linkage analysis, in which paternal and maternal transmissions of marker alleles are equally weighted. Several methods have been proposed for taking genomic imprinting into account in the model-free linkage analysis of binary traits. However, none of these methods are suitable for the formal identification and quantification of genomic imprinting in the presence of linkage. In addition, the available methods are designed for use with pure sib-pairs, requiring artificial decomposition in cases of larger sibships, leading to a loss of power. We propose here the maximum likelihood binomial method adaptive for imprinting (MLB-I), which is a unified analytic framework giving rise to specific tests in sibships of any size for (i) linkage adaptive to imprinting, (ii) genomic imprinting in the presence of linkage, and (iii) partial versus complete genomic imprinting. In addition, we propose an original measure for quantifying genomic imprinting. We have derived and validated the distribution of the three tests under their respective null hypotheses for various genetic models, and have assessed the power of these tests in simulations. This method can readily be applied to genome-wide scanning, as illustrated here for leprosy sibships. Our approach provides a novel tool for dissecting genomic imprinting in model-free linkage analysis, and will be of considerable value for identifying and evaluating the contribution of imprinted genes to complex diseases.

Aspects of genetic susceptibility to human infectious diseases.

Host genetic factors play a major role in determining differential susceptibility to major infectious diseases of humans, such as malaria, HIV/AIDS, tuberculosis, and invasive pneumococcal disease. Progress in identifying the relevant genetic loci has come from a variety of approaches. Most convincing associations have been identified by case-control studies assessing biologically plausible candidate genes. All six of the genes that have a major effect on infectious disease susceptibility in humans have been identified in this way. However, recently genome-wide linkage analysis of affected sibling pairs has identified susceptibility loci for chronic infections such as leprosy and chronic hepatitis B virus persistence. Other approaches used successfully have included assessment in humans of the homologues of susceptibility genes mapped and identified in murine models. However, the great majority of susceptibility loci remain to be identified and the advent of large-scale genome-wide association scans offers a new approach to defining many of these.

Genome-wide scan for loci influencing quantitative immune response traits in the Belém family study: comparison of methods and summary of results.

Here we report the results from a genome-wide linkage scan to identify genes and chromosomal regions that influence quantitative immune response traits, using multi-case leprosy and tuberculosis families from north-eastern Brazil. Total plasma IgE, antigen-specific IgG to Mycobacterium leprae soluble antigen (MLSA), M. tuberculosis soluble antigen (MTSA) and M. tuberculosis purified protein derivative (PPD), and antigen-specific lymphocyte proliferation (stimulation index or SI) and interferon-gamma (IFN-gamma) release to MLSA and PPD, were measured in 16 tuberculosis (184 individuals) and 21 leprosy (177 individuals) families. The individuals were genotyped at 382 autosomal microsatellite markers across the genome. The adjusted immune-response phenotypes were analysed using a variety of variance components and regression-based methods. These analyses highlighted a number of practical issues and problems with regard to implementation of the methods and, interestingly, differences were observed between several standard statistical and genetic analysis packages used. From this we determined that, for this set of traits in these pedigrees, significant p values for linkage using variance components analysis, supported by significance using the Visscher-Hopper modification of the Haseman-Elston method, provided the most compelling evidence for linkage. Using these criteria, linkage (5.8 x 10(-5) < p < 0.008) was seen for: total plasma IgE on chromosome 2; IgG to MLSA on chromosomes 8, 17 and 21; IgG to PPD on chromosome 12; SI to PPD on chromosome 1; IFN-gamma to MLSA on chromosomes 6, 7, 10, 12 and 14; and IFN-gamma to PPD on chromosomes 1, 16 and 19.

Linkage of leprosy susceptibility to Parkinson's disease genes.

In early 2003, an international team of scientists conducted a genome scan in Vietnamese multiplex leprosy families and found that susceptibility to leprosy was significantly linked to region q25 on the long arm of chromosome 6. Further confirmation of the chromosome 6 locus was provided by high-resolution linkage mapping in simplex leprosy families. Now, in a continuation of these findings, the team has pinpointed the chromosome 6 susceptibility locus to the 5' regulatory promoter region shared by both the Parkinson's disease gene PARK2 and its co-regulated gene PACRG. The surprising discovery has important implications for the understanding of leprosy pathogenesis and for the strategy of genetic analysis of infectious diseases.

Linkage analysis of susceptibility to leprosy type using an IBD regression method.

Leprosy is a chronic disease caused by infection with Mycobacterium leprae, which is manifested across a wide clinical spectrum. There is evidence that susceptibility both to leprosy per se and to the clinical type of leprosy is influenced by host genetic factors. This paper describes the application of an identity by descent regression search for genetic determinants of leprosy type among families from Karonga District, Northern Malawi. Suggestive evidence was found for linkage to leprosy type on chr 21q22 (P<0.001). The methodological implications of the approach and the findings are discussed.